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ABSTRACT 

Sympathetic eruptions on the Sun have been observed for several decades, but the mechanisms by which one 
eruption can trigger another one remain poorly understood. We present a 3D MHD simulation that suggests 
two possible magnetic trigger mechanisms for sympathetic eruptions. We consider a configuration that con- 
tains two coronal flux ropes located within a pseudo-streamer and one rope located next to it. A sequence of 
eruptions is initiated by triggering the eruption of the flux rope next to the streamer. The expansion of the rope 
leads to two consecutive reconnection events, each of which triggers the eruption of a flux rope by removing 
a sufficient amount of overlying flux. The simulation qualitatively reproduces important aspects of the global 
sympathetic event on 2010 August 1 and provides a scenario for so-called twin filament eruptions. The sug- 
gested mechanisms are applicable also for sympathetic eruptions occurring in other magnetic configurations. 
Subject headings: Sun: corona — Sun: coronal mass ejections (CMEs) — Sun: flares — Sun: filaments, 
prominences — Sun: magnetic topology — Methods: numerical 



1. INTRODUCTION 

Solar eruptions are observed as filament (or prominence) 
eruptions, flares, and coronal mass ejections (CMEs). It is 
now well established that these three phenomena are differ- 
ent observational manifestations of a single eruption, which 
is caused by the destabilization of a localized volume of the 
coronal magnetic field. The detailed mechanisms that trigger 
and drive eruptions are still under debate, and a larg e num- 
ber of theoretical models have been developed (e.g. 
[20T0] >. 

Virtually all existing models consider single eruptions. The 
Sun, however, also produces sympathetic eruptions, which oc- 
cur within a relatively short period of ti me - either in one, typ- 
ically complex, active region (e.g., Liu et al.|2009| > or in dif- 
ferent source regions, which occasionally cover a full hemi- 
sphere (so-called "global eruptions"; Zhukov & Veselovsky 



Forbes 



2007 ). It has been debated whether the close temporal cor- 



relation between sympathetic eruptions is purely coinciden- 
tal, or whether t hey are causally linked (e.g., |Bie seeker &| 
Thompson 2000 ). Both statistical i nvestigations (e.g., Moon 
et al.|2 002 ; Wheatland & Craig 2006 ) and detailed case stud- 
ies (e.g., |Wang et al.|2001| l indicate that physical connections 
between them exist 1 . 

The exact nature of these connections has yet to be estab- 
lished. They have been attributed, for instance, to convective 
motions or destabilization by large-scale waves (e.g., Ram- 



|sey & Smith][p66| IBumba & Klvana|[T993] ). At present, it 
seems most likely that the mechanisms by which one erup- 
tion can trigger another one act in the corona and are of a 
magnetic nature. Perturbations traveling alo ng field fines that 
connect source regions of eruptions (e.g., Jiang et aL]|2008 



and changes in the background field du e to reconnection (e.g., 
|Ding et al.|2006|[Zuccarel lo et al. 2009 ) have been considered. 
In an analysis of a global sympathetic event (see Section|2]i, 

1 We do not distinguish here between sympathetic flares and sympathetic 
CMEs, since both are part of the same eruption process. 



Schrijver & Title (2011 1 found evidence for connections be- 
tween all involved source regions via structural featu res like 
separators and quasi-separatri x layers (QSLs; [Priest & Forbes| 
1992| Demoulin et al. |l996| l; suggesting the importance of 
the structural properties of the large-scale coronal field in the 
genesis of sympathetic eruptions. 

A magnetic configuration that appears to be prone to 
producing sympathetic erupti ons is a unipolar streamer or 
pseud o-streamer (PS; e.g., Hundhausen] 1 1 972] Wang et al. 
|2007| l. A PS is morphologically similar to a helmet streamer, 
but divides open fields of like polarity and contains an even 
number (typically two) of closed flux lobes below its cusp. 



PSs are quite common in the corona (e.g., |Eselevich et al. 
1999| Riley & Luhmann]|2011|l and occasionally harbor two 



filaments. It seems that if one of these erupts, the other 
one fo llows shortly thereafter ( so-called "twin filament erup- 
tions"; |KnaseMO^Wiil20T0|. 

Here we present a numerical simulation that suggests two 
possible magnetic trigger mechanisms for sympathetic erup- 
tions. It was inspired by the global sympathetic event on 2010 
August 1, which involved a twin filament eruption in a PS. 

2. THE SYMPATHETIC ERUPTIONS ON 2010 AUGUST 1 

A detailed account of the individual eruptio ns that oc- 
curred in this global event can be found in|Schrijver & Ti- 



lde] ( |201 l| l. Here we focus on a subset of three consecu- 
tive filament eruptions, all of which evolved into a sepa- 
rate CME. FiguresfTL b, and c sh ow, respectively, the e rup- 
tions as seen by STEREO/EUVI ( [Howard et al]|2"0"08l >, the 
pre-emptive filamen ts, and a synoptic m agnetogram obtained 
from SoHO/MDI ( |Scherrer et al][T995] > data. The large fil- 
aments 2 and 3 were located along the inversion lines divid- 
ing an elongated positive polarity and two bracketing negative 
polarities; the small filament 1 was located at the edge of the 
southern negative polarity . A potential field so urce surface 
extrapolation (PFSS; e.g., [Schatten et al.||1969j l for Carring- 
ton rotation 2099 reveals that filaments 2 and 3 were lo cated 
in the lobes of a PS (Figure[T]l; see also |Panasenco & Velh| 
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FIG. 1. — (a:) STEREO-Ahead/EUVI 304 A observations of three subsequent prominence eruptions (marked by their order of eruption) on 2010 August 1, 
shown at 02:56, 09:16, and 22:06 UT (from left to right), (b): Big Bear Observatory Ha observation on 2010 July 30, showing the corresponding pre-emptive 
filaments, (c): Filament contours (drawn by eye) overlaid on a synoptic magnetogram for Carrington rotation 2099, with red (blue) showing positive (negative) 
radial fields, (d): Magnetic field lines from a corresponding PFSS extrapolation, revealing a pseudo-streamer. Green lines outline the lobes in which filaments 2 
and 3 were located, pink lines show adjacent coronal holes, (e): Coronal distribution of Q (grayscale) and photospheric distribution of slog Q, where red (blue) 
outlines positive (negative) magnetic fluxes. 



[20T0l l. 

Figure[lJ shows a cut through the coronal distribution of the 
squashing factor Q ( Titov et al.|2002| > above filaments 2 and 3. 
The dark lines of high Q outline structural features and exhibit 
here a shape characteristic for a PS (compare with Figure^ 
below). Th e photospheric distribution shows slog Q ( |Titov| 



et al.|201 T) , depicting the footprints of (quasi-)separatrix sur- 
faces. Trie structural skeleton of a PS consists of two separa- 
trix surfaces, one vertical a nd one dome-like, which are both 
surrounded by a thin QSL (|Masson et al.|2009| l and intersect 
at a separator ( |Titov et al.|201 l| l. It has been demonstrated that 
current sheet formati on and reconnection occur preferably at 



such separators (e.g., Baum & Bratenahl] 1980 Lau & Finn 
[T990> . 

The presence of the PS above filaments 2 and 3 suggests 
that the CME associated with filament eruption 1 may have 
triggered the subsequent eruptions by destabilizing the PS, 
presumably by inducing reconnection at its separator. We now 
describe an MHD simulation that enabled us to test this sce- 
nario using an idealized model. 

3. NUMERICAL SIMULATION 



The basic simulation setup is as in Toro k et al.| ( |201 l|l, 
where two instances of the coronal flux rope model by |Titov| 
|& Demoulin ( |1999[ hereafter TD) were used to simulate the 
interaction of two flux ropes in a PS. Here we add a TD con- 
figuration on each side of the PS (Figure[2k). The new con- 
figuration on the left is used to model the CME associated 
with filament eruption 1, while the new one on the right is 
merely used to obtain a (line-)symmetric initial configuration, 
which facilitates the construction of a numerical equilibrium. 
It does not significantly participate in the dynamic evolution 
described below. The flux ropes FR1-3 are intended to model 
filaments 1-3. 



We integrate the zero /3 compressible ideal MHD equations, 
neglecting thermal pressure and gravity. The equations are 
normalized by the initial TD torus axis apex height, R — d 
(see TD), the maximum initial magnetic field strength and 
Alfven velocity, Bomax and v a o m ax, and derived quantities. 
The Alfven time is r a = (R — d)/v a o max . We use a nonuni- 
form cartesian grid of size [—25, 25] x [—25, 25] x [0, 50] with 
resolution ~ 0.04 in the flux rope area. The initial density dis- 
tribution is po(x) = |£?o (x)\ 3 / 2 , such that v a (x) decreases 
slowly with distance from th e flux concentrations. F or further 
numerical details we refer to Torok & Kliem (2003 ). 

The model parameters are chosen such that al l flux ropes 
are in itially stable with respect to the helical k ink (Toro k et al.| 
|20031 l and torus instabilities (TI; |Kliem & Torok|2006| ). The 
ropes are placed along the y direction, at x = ±1.5 and ± 
5.5, and have identical parameters (R — 2.75, a — 0.8, d — 
1.75, L = 0.5, q = 4.64; see TD). The signs of the sub- 
photospheric point charges, ±q, are set according to the signs 
of the polarities surrounding filaments 1-3 (FigurefT];). The 
half-distance between the charges, L, is such that the TI can 



be trig gered by a relatively weak perturbation ( Schrijver et al. 
2008 ). To obtain a numerically stable initial configuration that 



contains (semi-)open field above the PS lobes, the two charges 
associated with each flux rope are adjusted to —0.55 q/0.65 q 
(for FR1 and FR4) and to -0.34 g/0.24 q (for FR2 and FR3). 
The twist is chosen left-handed for all ropes to a ccount for the 
observed d extral chirality of filaments 2 and 3 ( Panasenco & 
|Velli|2010| ). 

We first relax the system for 85 r a and reset time to zero. 
Then we trigger the eruption of FR1 by imposing localized 
converging flows at the bottom plane (as in Toro k et al.|201 1[ ), 
which slowly drive the polarities surrounding FR1 toward the 
local inversion line, yielding a quasi-static expansion of the 
rope's ambient field. The flows are imposed for 25 r a (in- 
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FIG. 2. — Snapshots from the simulation, showing magnetic field lines with fixed footpoints and the normal component of the magnetic field at the bottom 
plane, where red (blue) depicts positive (negative) fields. Orange lines belong to the flux ropes, green ones to the initial pseudo-streamer lobes, and pink ones 
to initially closed or (semi-)open overlying flux. Panel (a) shows the configuration after initial relaxation and (b-d) show the successive flux rope eruptions and 
ambient field evolution at t = 85, 126, and 181 T a , respectively. "Already erupted" flux ropes are omitted for clarity. (An animation of this figure is available in 
the online journal.) 



eluding phases of linear increase [decrease] to a maximum 
velocity of 0.02 n o( | max [to zero], each lasting 5 r Q ). 

Though we solve the ideal MHD equations, extra diffusion 
is introduced by numerical differencing (as in every MHD 
code that models solar magnetic fields). This numerical diffu- 
sion is localized in regions where the current density is largest, 
and leads to reconnection of magnetic field lines. Although it 
is much larger than the diffusion expected on the Sun, experi- 
ence has shown that simulations produce solutions with phys- 
ically expected behavior, as long as the numerical diffusion 
is sufficiently small. We therefore expect that our simulation 
indicates the true evolution of the system, but that the recon- 
nection rates might differ from those present on the Sun. 

4. RESULTS 

Figures[2] and [3] summarize the main dynamics and recon- 
nection occurring in the simulation. Figure[3^ shows the ini- 
tial configuration and Figures^ and show the system af- 
ter relaxation, during which weak current layers form at the 
PS separatrix surfaces, but no noticeable reconnection occurs. 
Note the correspondence between the current layer pattern 
and the Q-distribution shown in Figureflt. 

As the converging flows are applied, FR1 starts to rise 
slowly, in response to the quasi-static expansion of its ambient 
field. In contrast to other simulations, where such flows have 
been used to crea te a flux rope from a sheared arcade (e.g., 
Amari et al.|2000| l, they do not lead here to noticeable recon- 
nection. The slow rise lasts until the rope reaches the critical 
height for TI onset at t ~ 40 r a , after which it rapidly accel- 
erates upward driven by the instability (fTorok & Kliem|2007 



2Q1Q| >. FR1 attains a maximum velocity of w 0.45i? a o ma x 
at t « 90 T a before it slowly decelerates. FigurelZb shows 
the system in the course of this eruption. The rise orthe rope 
is slig htly inclined, due to the asymmetry of its ambient field 
(e.g., Filippov et al. 2001). The rope rotates counterclock- 
wise about its rise direction (as seen from a bove), due to the 
conversion of its twist into writhe (e.g., Gre en et al.|2007[ ). 

The expansion of FRl's ambient field compresses the field 
between FR1 and the PS, particularly at larger heights where 
it is weak (see online animations). As a result, a tilted arc- 
shaped current layer forms around the PS separator (Fig- 
ures^ and[4}. Further compression by the eruption steepens 
the current densities until reconnection (Rl) between the open 
flux to the left of the PS and the closed flux in the right PS lobe 
sets in. The lobe flux then starts to open up, while the open 
flux starts to close down above the left PS lobe (Figures^ and 
[3J;). This successively decreases (increases) the magnetic ten- 
sion above FR2 (FR3), so that FR2 rises slowly, while FR3 is 
slowly pushed downward. At t w 95 r a FR2 reaches the criti- 
cal height for TI onset and erupts, attaining a maximum veloc- 
ity of ~ 0.60 WaOmax at i « 120 r a . Figure[2l shows that FR2 
also rises non-radially, but rotates less thanFRl. The appar- 
ently smaller rotation of FR2 is due to the faster decay of its 
overlying field with height, which leads to a distributio n of the 
total rotatio n over a larger height range than for FR1 ( ]Tor6k| 
et al.||2010) l. By the time shown, FR1 has fully erupted, an 



|Fan & Gibsoii|[2007| |Schrijver etaLl|2008[ |Aulanier et al. 



elongated vertical current layer has formed in its wake (Fig- 
ure[3]l), and reconnection therein has produced cusp-shaped 
fiela lines below it. As FR2 erupts, it rapidly pushes the arc- 
shaped current layer to large heights (Figure[3jl). While Rl 
still commences for some time, it does not play anymore a 
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FIG. 3 . — Illustration of the two reconnection phases that trigger the successive flux rope eruptions in the pseudo- streamer. Field lines are colored as in Figure|2] 
The transparent inverted grayscale in the central plane, {y = 0}, shows the logarithmic distribution of |j|/|B|, where j is the electric current density, outlining 
flux rope currents and thin current layers. Fainter field line segments are located behind the transparent layer. Panel (a) shows the initial configuration, (b) the 
system after relaxation, and (c-f) show the dynamic evolution, at t = 85, 126, 142, and 158 r a , respectively. Panels (a,b) and (c-f) use a different scaling of 
|j|/|B|, respectively. Panel (c) shows reconnection Rl, which triggers the eruption of FR2, and (e,f) show reconnection R2, which triggers the eruption of FR3. 
Panel (d) shows a state between the two reconnection phases. (An animation of this figure is available in the online journal.) 



significant role for the following evolution. 

A vertical current layer also forms below FR2. The subse- 
quent reconnection (R2) initially involves the very same flux 
systems that took part in Rl . The flux previously closed down 
by Rl opens up again, and the flux previously opened up by 
Rl - and by the expansion of FR2 - closes down to form cusp- 
shaped field lines below the current layer (Figure[3£). After 
these fluxes are exhausted, R2 continues, now involving the 
left PS lobe and the open flux to the right of the PS. While the 
former opens up, the latter closes down as part of the growing 
cusp (Figure[3t). Thus, R2 continuously removes closed flux 
above FR3. As before, this progressive weakening of mag- 
netic tension leads to a slow rise of the rope, followed by its 
eruption (Figures^ andpjl). The rapid acceleration of FR3 
by the TI starts at t « 138T a , yielding a maximum velocity 
of w 0.35 VaOmax at t « 175 T a . The rope shows a significant 
rotation and an inclined rise which is now mainly directed to- 
ward the positive x direction. 

5. DISCUSSION 

The eruptions of FR2 and FR3 are initiated by the re- 
moval of a sufficient amount of stabilizing flux above the flux 
ropes via reconnection. Rl is similar to quadrup olar "break- 
through" or "break out" reconnection ( |Syrovatskii| 1 982| |Anti-| 
ochos et al.| 1999| >. Here it is driven by a nearby CME rather 
than by an expanding arcade and, in contrast to the breakout 
model, a flux rope is present prior to eruption. R2, on the other 
hand, corresponds to standard flare reconnection in the wake 
of a CME. Here it removes flux from the adjacent PS lobe, 
thereby triggering the eruption of FR3. A similar mechanism 
for the initiation of a s econd eruption in a PS was suggested 
by Cheng et al. ( 2005 i, who, however, attributed it to recon- 
nection inflows rather than to flux removal. We emphasize 
that Rl and R2 merely trigger the eruptions, which are driven 



by the TI and su pported by the associated flare reconnection 
(e.g.,[Vrsnak 2008 ). Thus, in the system studied here, both PS 
eruptions require the presence of a pre-emptive flux rope. We 
further note that the reconnections do not have to commence 
for the whole time period until the TI sets in. It is sufficient 
if they remove enough flux for the subsequently slowly rising 
flux ropes to reach the critical height for TI onset. 

Rl is driven by a perturbation of limited duration - the lat- 
eral expansion of a nearby CME - and is slow since it involves 
only weak fields, around a separator at a significant height in 
the corona. Therefore, its success in triggering an eruption 
depends on parameters like the distance of the CME from the 
PS and the amount of pre-emptive flux within the PS lobes. 
Indeed, if we sufficiently increase these parameters in the sim- 
ulation we find that Rl still commences, but does not last long 
enough to trigger an eruption. In contrast, R2 is driven by the 
rise of FR2 and involves strong fields. It is ther efore faster 
and more effic ient, which supports the finding by |Panasenco| 



& Velli (2010 1 that an eruption in one lobe of a PS is often 
followed by an eruption in the neighboring lobe. 

Figure |2]shows that the simulation correctly reproduces the 
order of the eruptions shown in Figure[T^ and yields a good 
match of their inclinations and rotations. Assuming that the 
first eruption indeed triggered the subsequent ones, it is sur- 
prising that the filament located further away from it went 
off first. While filament 2 may simply have been closer to 
its stability limit than filament 3 (as indicated by its larger 
height; see Figure[T^), the simulation provides an alternative 
explanation: the perturbation of the separator yields an ori- 
entation of the current layer that leads to a removal of closed 
flux only in the right PS lobe (Figure|3};), thus enforcing the 
eruption of FR2. Hence, although we could not find obser- 
vational signatures of Rl (presumably because the involved 
fields were too weak), the observed eruption sequence sup- 
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FIG. 4. — Volume renderingof |j|/|B| in the pseudo-streamer area at the 
same time as in Figures|2j) andFJJ:, outlining the tilted arc-shaped current layer 
that forms around the separator. 

ports its occurrence. The time intervals between the simu- 
lated eruptions exhibit a ratio different to the observed ones. 
Matching the observed ratio requires a search for the appro- 
priate model parameters and a more realistic modeling of re- 
connection, which are beyond the scope of this work. 

FR2 reaches a velocity about 35 per cent larger than FR1, 
which is in line with |Liu| ( |2007| l and [Fainshtei n & Ivanov| 
( |2010| ), who found that CMEs associated with PSs are, on 
average, faster than those associated with helmet streamers. 
|Liu| ( |2007|l suggested that this difference is due to the typi- 
cally smaller amount of closed flux the former have to over- 
come. Indeed, FR1 has to pass through flux that is closed at 
all heights above it, while FR2 faces much less closed flux, 
a significant fraction of which is, moreover, removed by Rl. 
FR3 remains significantly slower than FR2, most likely be- 
cause it encounters more closed flux at eruption onset, and 
only partially opened flux later on (Figures[3£,f). 

6. CONCLUSIONS 

We present an MHD simulation of two successive flux rope 
eruptions in a PS, and we demonstrate how they can be trig- 
gered by a preceding nearby eruption. The simulation sug- 
gests a mechanism for twin filament eruptions and provides 



a scenario for a subset of the sympathetic eruptions on 2010 
August 1 . More realistic initial configurations and a more so- 
phisticated treatment of reconnection are needed for a quanti- 
tative comparison with observations. 

Our results support the conjecture that the trigger mecha- 
nisms of sympathetic eruptions can be related to the structural 
properties of the large-scale coronal field. However, while 
structural features are present in our model configuration, they 
do no connect the source region of the first eruption with the 
source regions of the subsequent ones. Moreover, the mere 
presence of such features in a source region is not a suffi- 
cient criterion for the occurrence of a sympathetic event, even 
if reconnection at structural features is triggered by a distant 
eruption. The conditions in the source region must be such 
that the resulting perturbation forces the region to cross the 
stability boundary. 

The two trigger mechanisms presented here are indepen- 
dent and applicable also to other magnetic configurations. 
Triggering a sympathetic eruption by Rl requires the presence 
of a separator (or null point) above closed flux that stabilizes a 
pre-emptive flux rope, which can be realized, in t he simplest 
case, in a so-called "fan-spine" confi guration (e.g., Antiochos 
[T998l|Pariat et al.|2009l|Torok et al.|2009) . Triggering a sym- 
pathetic eruption by R2 requires the presence of an adjacent 
closed flux system overlying a flux rope, which can exist, for 
example, in quadrupolar configurations. 
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